format long;
clear;
CX=[1e-10;2e-10;3.3e-10;4.7e-10;5.6e-10;6.8e-10;8.2e-10;10e-10;
1e-9;2e-9;3.3e-9;4.7e-9;5.6e-9;6.8e-9;8.2e-9;10e-9;
20e-9;33e-9;47e-9;56e-9;68e-9;82e-9;100e-9;
200e-9;330e-9;470e-9;560e-9;680e-9;820e-9;1000e-9];
[h,v]=size(CX);
k1 = 1;
for i=1:h
  for j=1:h
    C1 = CX(i);
    C2 = CX(j);
    a=C2-4*C1;
    b=-3e-4;
    c=1e-8/(C1);
    dt = b*b-4*a*c;
    if(dt > 0)
      R1=(-b - sqrt(dt))/(2*a);
      R2 = 1e-8/(R1*C1*C2);
      if((R1>500) && (R1 < 15e3) && (R2>500) && (R2 < 15e3))
        CR(k1, 1)=C1;
        CR(k1, 2)=C2;
        CR(k1, 3)=R1;
        CR(k1, 4)=R2;
        k1 = k1+1;
      end
    end
  end
end


